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A simple model of spacetime foam, made by N wormholes in a semiclassical 

(J., • approximation, is taken under examination. We show that the Bekenstein- 

T I^ • Hawking entropy is here "quantized" in agreement with the heuristic calcula- 



X 



tion of Bekenstein. 

I. INTRODUCTION 

In the early seventies J. Bekenstein |L| and S. Hawking ||, basing themselves upon con- 
siderations originated from quantum field theory, changed drastically the concept of a black 
hole. Beginning with the simple observation that the area of the horizon of the Schwarzschild 
black hole is a quadratic function of the mass M, they considered an infinitesimal increment 
(in natural units G — h — c — 1) 

dM = KdA hor k = — — . 



This formally resembles the First Law of thermodynamics dll = TdS^ This analogy of 
the black hole mechanics with thermodynamics appears reinforced by Hawking's theorem 
on black hole mechanics || asserting that the horizon area of an isolated black hole never 
decreases in any transformation. Thus, if two black holes of area A\ and A 2 fuse to form 
a black hole of area Ai +2 , then the Theorem asserts that A 1+2 > A\ + A 2 . On the basis of 
these observation and results, Bekenstein made the proposal that a black hole does have an 
entropy proportional to the area of its horizon 

Sbh = const x A hor . (1) 

In particular, in natural units one finds that the proportionality constant is set to 1/4G = 
1/4/p, so that the entropy becomes 

For a Schwarzschild black hole, for example, one finds the value 



S = — !> 1- = AttM 2 G = AttM 2 E 

AG p 

In conventional units instead, we find for a generic horizon area that 



(3) 



S bh = ^-^r In 2k B c 3 A hor . (4) 



x The best way to see this similarity is with a Kerr-Newman black hole, where the incremental 
formula is 

dM = KdA hor + &dQ + O • dt 

and the analogy with the First Law 

dU = TdS + PdV 

is complete. 



This formula is identical (except for the factor of In 2 which one may think of as a choice of 
units of entropy) to the one proposed by Hawking || based on consistency with the rate of 
black hole radiation. The appearance of h is "... a reflection of the fact that the entropy is 
... a count of states of the system" ]]]]. Following Bekenstein's proposal on the quantization 
of the area for nonextremal black holes we have 

a n = alp (n + rf) r\ > — 1 n = 1, 2, . . . 

Many attempts to recover the area spectrum have been done, see Refs. MM for a review. Note 
that the appearance of a discrete spectrum is not so trivial. Indeed there are other theories, 
based on spherically symmetric metrics in a mini-superspace approach, whose mass spectrum 
is continuous. Recently a model made by N coherent wormholes, based on Wheeler's ideas 



of a foamy spacetime ||, has been considered |12| . In that paper, hereafter referred as I, 
we have computed the energy density of gravitational fluctuations reproducing the same 
behavior conjectured by Wheeler during the sixties on dimensional grounds. In this paper 
we wish to apply the ideas of I to a generic area. The result is a quantization process whose 
quanta can be identified with wormholes of Planckian size. Implications on the black hole 
entropy are taken under consideration. The rest of the paper is structured as follows, in 
section [TT] we briefly recall the results reported in Ref. ||12| , in section [TlT] we compute the 
area spectrum. We summarize and conclude in section [TV] . 

II. SPACETIME FOAM: THE MODEL 

In the one-wormhole approximation we have used an eternal black hole, to describe a 
complete manifold A4, composed of two wedges M.+ and M.- located in the right and left 
sectors of a Kruskal diagram. The spatial slices £ represent Einstein-Rosen bridges with 
wormhole topology S 2 x R 1 . Also the hypersurface £ is divided in two parts £+ and £_ by 
a bifurcation two-surface Sq. We begin with the line element 

ds * = _at2 ( r ) d ? + _^!_ + r 2 (d6 2 + sin 2 6d<p 2 ) (5) 



and we consider the physical Hamiltonian defined on £ 

H P = H-H = ^—J d 3 x (NH + N{W) 

+ ^J d 2 x^(k-k°)-^J d 2 x^(k-k°), (6) 

where It = G. The volume term contains two constraints 

H = G^H- (|) - (f ) J*» = ^ ( ^ 

W« = -2vrg = 

both satisfied by the Schwarzschild and Flat metric respectively. The supermetric is Gijki = 
2 (dikdji + 9u9jk — 9ij9ki) and i?® denotes the scalar curvature of the surface E. By using 



the expression of the trace 



k 



^(Vhn») , (8) 



Vh 

i 
with the normal to the boundaries defined continuously along S as n M = (h yy ) 2 S£. The 

value of k depends on the function r, y , where we have assumed that the function r, y is 

positive for S + and negative for SL. We obtain at either boundary that 

— 1r 
k = -=^L. (9) 

r 

The trace associated with the subtraction term is taken to be k° = — 2/r for B + and k° = 2/r 
for B_. Then the quasilocal energy with subtraction terms included is 

£ qua siiocai = ? (E+ ~ E.) = ll [(r [1 - \r, y \}) y=y+ - (r [1 - \r, v \}) y=y _] . (10) 

Note that the total quasilocal energy is zero for boundary conditions symmetric with respect 
to the bifurcation surface Sq and this is the necessary condition to obtain instability with 
respect to the flat space. In this sector satisfy the constraint equations (|7|). Here we consider 
perturbations at S of the type 

9ij = 9ij +hij, (11) 



where gij is the spatial part of the Schwarzschild and Flat background in a WKB approxi- 
mation. In this framework we have computed the quantity 



AE(M) 



* 



H 



Schw. 



-H 



Flat 



* 



{*\H, 



quasilocal 



*) 



(12) 



by means of a variational approach, where the WKB functionals are substituted with trial 
wave functionals. This quantity is the natural extension to the volume term of the subtrac- 
tion procedure for boundary terms and it is interpreted as the Casimir energy related to 
vacuum fluctuations. By restricting our attention to the graviton sector of the Hamiltonian 
approximated to second order, hereafter referred as H\ 2 , we define 



E 



|2 



#" 



H\ 



y- 



(tf-L|tf-L) ' 



where 



m x = V 



Kl 



jVexp 



^ 



(ig-g)K' l {g-g)f 

\ I x,y 



After having functionally integrated H\ 2 , we get 



4/Hjs 



H l2 = — dtxy/gCT"* K~^ (x,x) ijkl + (A 2 ) a j K ± {x ) x) l 



i iakl 



(13) 



The propagator K x (x, x) iakl comes from a functional integration and it can be represented 

as 



^(^i/W=E 



JV 



Ha (X) h kl (it) 
2X N (p) 



(14) 



where hj~ a (a?) are the eigenfunctions of 



'A 2 r := -A<J- + 2m. 



(15) 



This is the Lichnerowicz operator projected on S acting on traceless transverse quantum fluc- 
tuations and Ajv (p) are infinite variational parameters. A is the curved Laplacian (Laplace- 
Beltrami operator) on a Schwarzschild background and W is the mixed Ricci tensor whose 
components are: 



nn ,. (-2MG MG MG) 

^ = ^|— ^,-r,-rj. (16) 

The minimization with respect to A and the introduction of a high energy cutoff A give to 
the Eq. flT2p the following form 

A^/,^ y /3MG\ 2 , /r 3 A 2 \ 

A£ (M) ~ — 5— In — — , (17) 

V ; 32tt 2 \ rg / \SMGJ ' V ; 

where V is the volume of the system and ro is related to the minimum radius compatible with 

the wormhole throat. We know that the classical minimum is achieved when r$ = 2MG. 

However, it is likely that quantum processes come into play at short distances, where the 

wormhole throat is defined, introducing a quantum radius r > IMG. We now compute 

the minimum of AE(M), after having rescaled the variable M to a scale variable x = 

3MG/{r 3 A 2 ). Thus 

AE (M) -> AE (x, A) = -A 4 i 2 lnx 

v ; v ' ; 327T 2 

We obtain two values for x: X\ = 0, i.e. flat space and x<i = e~z. At the minimum 

AE< ^ = -^T- (18) 

Nevertheless, there exists another part of the spectrum which has to be considered: the 
discrete spectrum containing one mode. This gives the energy an imaginary contribution, 



namely we are discovered an unstable mode OJTBI]. Let us briefly recall, how this appears. 
The eigenvalue equation 

(A 2 ) a t h aj = ah tJ (19) 

can be studied with the Regge- Wheeler method. The perturbations can be divided in odd 
and even components. The appearance of the unstable mode is governed by the gravitational 
field component h e ^ n . Explicitly 

-E 2 H (r) 

2MG\d 2 H(r) (It - 3MG\ dH (r) AMG , , 



where 



H (r) 1 



ht? n (r,#,<f>) 



and E 2 > 0. Eq. (EUf) can be transformed into 



2m ^ - ] 
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Yoo^A) 



(21) 
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/dM^y 



My)' 



Mv) 



I dyh 2 (y) 



(22) 



where /i is the eigenvalue, y is the proper distance from the throat in dimensionless form. If 
we choose h(X,y) = exp {—Xy) as a trial function we numerically obtain /i = —.701626. In 
terms of the energy square we have 



E 2 = -.175 41/ (MGf 



(23) 



to be compared with the value E 2 = — . 19/ (MG) of Ref. JTT||. Nevertheless, when we 
compute the eigenvalue as a function of the distance y, we discover that in the limit y — > 0, 



fi= ^(A) = A 2 -- + - 



-2 . y 
y + Yx 



(24) 



Its minimum is at A = (■^y) 3 and 



/i (T) = 1.2878^ + 



- 9 _ 2 3 



-y 



(25) 



It is evident that there exists a critical radius where fi turns from negative to positive. This 
critical value is located at p c — 1. 113 4 to be compared with the value p c = 1.445 obtained by 



B. Allen in ||17|| . What is the relation with the large number of wormholes? As mentioned 
in I, when the number of wormholes grows, to keep the coherency assumption valid, the 
space available for every single wormhole has to be reduced to avoid overlapping of the wave 
functions. If we fix the initial boundary at R±, then in presence of iV„, wormholes, it will 
be reduced to R±/N w . This means that boundary conditions are not fixed at infinity, but 
at a certain finite radius and the ADM mass term is substituted by the quasilocal energy 
expression under the condition of having symmetry with respect to each bifurcation surface. 



The effect on the unstable mode is clear: as N w grows, the boundary radius reduces more 
and more until it will reach the critical value p c below which no negative mode will appear 
corresponding to a critical wormholes number N Wc . To this purpose, suppose to consider 
N w wormholes and assume that there exists a covering of E such that E = U Ej, with 

8=1 

Ej fl Ej = when i ^ j. Each Ej has the topology S 2 x R 1 with boundaries <9Ef with 
respect to each bifurcation surface. On each surface Ej, quasilocal energy gives 

Thus if we apply the same procedure of the single case on each wormhole, we obtain 

E 



Of 9 /■ 

Ei quasilocal = « / (Pxy/v (k ~ fc°) ~ J^ / C^Xy/ff (k - fc°) . (26) 



\ quasilocal = ? (^+ " £*-) = § (r [1 - \r, y \}) y=y%+ ~ jj {v [1 - \r jy |])^_ . (27) 

Note that the total quasilocal energy is zero for boundary conditions symmetric with respect 
to each bifurcation surface Sq^. We are interested in a large number of wormholes, each of 
them contributing with a term of the type (|T2"D. If the wormholes number is N w , we obtain 
(semi classically, i.e., without self-interactions)0 

ffffi = H 1 + H 2 + . . . + H Nw r (28) 

Thus the total energy for the collection is 

E\2 = N W H\ 2 . 
The same happens for the trial wave functional which is the product of N w t.w.f.. Thus 



2 Note that at this approximation level, we are in the same situation of a large collection of N 
harmonic oscillators whose hamiltonian is 






TVexp 



(g-g)K- l (g-g) 



.i.y 



where we have rescaled the fluctuations h — g — g in such a way to absorb N w /ti. The 
propagator K 1 - (x, x) iakl is the same one for the one wormhole case. Thus, repeating the 
same steps of the single wormhole, but in the case of N w wormholes, one gets 

V 



AE Nw (x,A.)~N 2 -— AVlnx 



32tt 2 



(29) 



where we have defined the usual scale variable x = 3MG/ (rgA 2 ). Then at one loop the 
cooperative effects of wormholes behave as one macroscopic single field multiplied by A" 2 , 
but without the unstable mode. At the minimum, x 

V A 4 



_ i 
e 2 



AE(x) 



-N z 



64vr 2 



(30) 



III. AREA SPECTRUM AND ENTROPY 

A very important application of the model presented in the previous section is the area 
quantization. The area is measured by the quantity 



A(S 



So 



d 2 x\fa. 



(31) 



a is the two-dimensional determinant coming from the induced metric a^ on the boundary 
So- We would like to evaluate the mean value of the area 



A (S ) 



*j 



.4 



*j 



*; 



Isa d 2 x^/a 



*j 



computed on 



{*f\Vf) 



\\& F ) = \J/^- <g> ^ 



(Vf\Vf) 



(32) 



* 



N« 



(33) 



If we use the fact that a a h = a a b + 5cr a b, where a a b is such that J s d 2 x\fa = 4irr 2 and r is 
the radius of So, we obtain to the lowest level^ in the expansion of a a b that 



5 For lowest level in the expansion of <7 a &, we mean 



A (S ) 



ty F A ^ F 



(*f\Vf) 



A-nr 1 



(35) 



Suppose to consider the mean value of the area A computed on a given macroscopic fixed 

N 

radius R. On the basis of our foam model, we obtain A = \J Ai, with Ai fl Aj = when 

i=l 

i ^ j. Thus 



A? 



N 



A = AnR 2 = £ A = J2 4 ^ 



(36) 



i=\ i=\ 



When Ti — > l p , Ai —> A ip and 



A = NA,„ = NAttP. 



(37) 



Thus the macroscopic area is represented by iV microscopic areas of the Planckian size: 
in this sense we will claim that the area is quantized. One immediately observes that iV 
is such that iV > N Wc , where N Wc is the critical wormholes number above which we have 
the stability of our foam model. At this point we can apply the same reasoning of Refs. 
|18| , |19| . |2T| to arrive at the Bekenstein-Hawking relation between entropy and area 

.4 



S 



K 



(38) 



Note that there is a In 2 numerical factor missing to complete the equality, between this 



model and the models described in Refs. ||18lJl9| , |21| . This is principally due to the fact that 
there is no a degeneracy factor related to the statistics. Here the wormholes are the same. 
Of course the introduction of a degeneracy factor does not alter the main result and the 



fa^ h = exp Tr In y/a^, = exp Tr- In (a ab + <5cr a6 ) 



exp Tr— 



In a ab 1 + 



6a. 



ab 



Vab 



expTr— 



In a ab + In 1 + 



Sa„ 



&ab 



/Oab + O (5<J ab ) 



(34) 



10 



only effect will be that of dividing the spacetime covering into equivalence classes with a 
representative for each class. In our case, we deduce that the entropy is 

S = Ntt. (39) 

Moreover this seems also to agree with the conclusions of Bekenstein ( || and Refs. therein), 
apart the degeneracy factor that in our model seems to be related to the odd or even 
permutation of any single wormhole wavefunctionQ. We can use Eq. (J37f) to compute the 
entropy also for other geometries, for example, the de Sitter geometry. Since we know that 



for this metric the Bekenstein-Hawking relation (]38|) still holds, we write 



3tt A NAitll , s 

that isQ 

Thus the cosmological constant A is " quantized" in terms of l p . Note that when the worm- 
holes number N is quite "large" , A — ► 0. We could try to see what is the rate of change 
between an early universe value of the cosmological constant and the value that we observe. 
In inflationary models of the early universe is assumed to have undergone an early phase with 
a large effective A ~ (10 10 - lO n GeV) 2 for GUT era inflation, or A ~ (10 16 - 10 18 GeV) 2 for 



4 Remark. Entropy is "quantized" as a consequence of ([37|) and not as a direct application of the 
definition 5 = — yjp n lnp n . 
5 A relation relating A and G, via an integer N appeared also in Ref. [20|. Nevertheless in Ref. 



[20], N represents the number of scalar fields and the bound from above and below 

|2C7iVA/3-2| > V3 
comes into play, instead of the equality ( PH ) . 



11 



Planck era inflation. A subsequent phase transition would then produce a region of space- 
time with A < (10~ 42 GeV) , i.e. the space in which we now live. For GUT era inflation, we 
have (we are looking only at the order of magnitude) 

10 20 - 10 22 GeV 2 = -10 38 GeV 2 ^N= 10 16 - 10 18 , (42) 

while for Planck era inflation we have 

10 32 - 10 36 GeV 2 = -Uo 38 GeU 2 -> N = 10 6 - 10 2 , (43) 

to be compared with the value of (10~ 42 GeV) which gives a wormholes number of the order 
of 

KT 84 GW 2 = ll0 38 GW 2 - N = 10 122 . (44) 

This very huge number is obtained by averaging the area on a cosmological scale by means 
of Planck scale wormholes. Thus it seems quite reasonable that with the growing of the 
cosmological radius, we obtain a growing wormholes number covering the horizon area. 

IV. CONCLUSIONS 

In this paper we have applied the model presented in I to the entropy computation assum- 
ing the validity of the Bekenstein-Hawking relation. In this picture the area is "quantized'' in 
the sense that spacetime can be filled by a given integer number of disjoint non-interacting 
wormholes. This result is in agreement (apart a numerical factor) with the quantized area 
proposed heuristically by Bekenstein and also with the loop quantum gravity predictions of 
Refs. [ELRITJI], apart the degeneracy factor missing, principally due to the fact that we have 



an " 'ideal Boltzmann gas" of wormholes. A similar result appeared in Ref. [21]. Neverthe- 



less, in Ref. ||21|| , spacetime was assumed ab initio built up of cells of Planckian size, while 
it seems that, in order to have stability, spacetime needs to be covered by N wormholes of 
the Planckian size. Nevertheless, even if in this letter also the meaning of the cosmological 
constant seems to be related to a pure gravitational effect enforcing the idea that A is not 

12 



fundamental but is the effect of quantum fluctuations of the pure gravitational field, giving 
therefore strong indications of a foamy spacetime, we are not in presence of a model which 
can cure the well known problems of the still absent theory of quantum gravity. The main 
reason resides in a cutoff dependent model. Even if one can argue that a Planck length 
cutoff is quite natural, it is not clear how to compute such a value from first principles, even 
if it is likely that a relation between this formulation of quantum gravity based on a large 
number of coherent wormholes and the well accepted string theory could exist. 
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